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Abstract. We have studied two specific models of frustrated and disordered 
coupled Kuramoto oscillators, all driven with the same natural frequency, in the 
presence of random external pinning fields. Our models are structurally similar, 
but differ in their degree of bond frustration and in their finite size ground state 
properties (one has random ferro- and anti-ferromagnetic interactions; the other 
has random chiral interactions). We have calculated the equilibrium properties of 
both models in the thermodynamic limit using the replica method, with emphasis 
on the role played by symmetries of the pinning field distribution, leading to 
explicit predictions for observables, transitions, and phase diagrams. For absent 
pinning fields our two models are found to behave identically, but pinning 
fields (provided with appropriate statistical properties) break this symmetry. 
Simulation data lend satisfactory support to our theoretical predictions. 



1. Introduction 

The dynamics and the analysis of models of interacting oscillators have received an 
increasing amount of attention during the last few decades (l). These models are 
excellent for studying synchronization phenomena and other features related to the 
temporal activity of populations of microscopic elements with limit cycle behaviour. 
They have been studied in many fields, including physics, chemistry and biology 
One of the most successful models was proposed by Kuramoto ||, who assumed 
that each member of the population can be modeled as an oscillator moving in a 
globally attracting limit cycle of constant amplitude. The interactions between the 
oscillators are sufficiently weak to ensure that perturbations will not move them away 
from their individual limit cycles. Only one degree of freedom, the phase Oi, is then 
required to describe the dynamics of each oscillator i. The latter evolves according to 
2£0i — u>i + 52 ■ Kijf(9i — 9j). Here cjj is the natural frequency of oscillator i (drawn 
randomly from some distribution), is the coupling strength between oscillators 
i and j, and f(x) a non- linear function. Kuramoto initially considered mean field 
coupling strengths = K/N, where N is the total number of oscillators, and 
f(x) = sin(x). This particular model has been investigated intensively. Generally, 
for some critical (positive) value of the coupling K a phase transition occurs from a 
state where all the oscillators run incoherently to a state where a certain degree of 
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synchronization emerges spontaneously. The details of this transition depend on the 
the distribution of natural frequencies |[ ^, [7j and the nature of the interactions 



(viz. short range interactions [|8|, |9j, random site disorder 10 1, or even more complex 
types |ll]]). Kuramoto-type models have in addition been studied extensively in the 
presence of noise j|, Q |l3| . 

The equations of the Kuramoto model have also been derived in the context of 
neural networks [|l4| 15 and Josephson junctions jl6|, [I?], [l8|. There, however, one 
usually obtains an extra term Aij in the argument of the nonlinear function: 

jO, =co t + J2 K ijf( e i ~ °3 + A v) (!) 

3 

The simplest case, Aij — a for all was studied in p9j . Here the degree 

of spontaneous synchronization was found to be maximal when a = 0, decreasing 
monotonically to zero as a — ► ^. A similar, but more complicated case is the subject 
of |2(| . However, the problem is far from solved for the case of quenched bond disorder. 
Here only simulation studies ^2| have been published so far; these emphasize the 
crucial role played by the disorder in determining the long time properties of the 
model, as measured by correlation functions. The situation is even more complex 
when random pinning fields are added, leading to further competition. 

The simplest models with randomness in both interactions and external fields are 
those where all oscillators have the same natural frequency, u>i — lu; here one can switch 
to a new basis and transform away the u>i. Moreover, if the matrix {-Ay} fulfills certain 
symmetry criteria, then (Q) can be written as a gradient descent process, enabling 
(upon adding Gaussian white noise) an equilibrium statistical mechanical analysis. 
Low dimensional versions of the resulting system have also been studied; they are 
equivalent to frustrated XY models, used to describe Josephson junction arrays. Here 
the sum over interacting pairs in the Hamiltonian is restricted to neighbours which 
form a plaquette, with 9i denoting the orientation of the spin at site i, and Aij denoting 
the bond angle, such that the plaquette sum X)(ij)g P iaq^«i = 27r/, where / is a 
measure of the frustration. For / = 1/2 one has the so-called fully frustrated model. 
Both the critical behaviour and the symmetries of these low dimensional systems are 
known to depend highly discontinuously on / [ p3| . Random external fields trigger 
local fluctuations of /, giving rise to a complex phenomenology which is still subject 
of study. It would be interesting to investigate whether such features remain in mean- 
field models. To our knowledge, this has not yet been studied; the effect of random 
fields has so far only been analyzed in systems without quenched disorder [E3 , or for 



m- vector spin glass models 25, |26| [27| (but for local uni-axial an- isotropic fields of a 
different nature than the ones discussed here). 

Our paper is structured as follows. We first define two mean-field models 
with disorder in both bonds and pinning fields, and their appropriate macroscopic 
observables. The first model is of a conventional type, with individual pairs of 
oscillators trying to either fully synchronize or fully anti-synchronize, depending on the 
sign of their interaction. The second model is less conventional in that neuron pairs 
prefer phase differences Aij of either tt/2 or — tt/2. For both models the presence of 
random external pinning forces increases energetic conflicts and frustration further. In 
section three we solve both models in equilibrium, using the replica method; we make 
the replica-symmetric ansatz and calculate the conditions for replicon instabilities. In 
sections four and five we study the effects of global symmetries. Finally, in section 
six we present results (in the form of phase diagrams and the temperature and field 
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strength dependence of observables) for a number of specific choices for the pinning 
field distributions, and validate our predictions via numerical simulations. 

2. Model definitions 

We study systems of N coupled Kuramoto oscillators (or XY spins) with external 
pinning fields, described by Hamiltonians of the form 
2K 

H= — =23cos(0(-0 J +^-) -^X/ C0S ^ - ^) ( 2 ) 

i<3 i 

where Oi S [0, 2n\ denotes the phase of the i-th oscillator. The natural dynamics 
leading to a Boltzmann state with Hamiltonian (^|) is the Langevin equation 4^9i — 
—dH/dOi + £i(t), with Gaussian random forces £j(t) which obey (£i(t)) = and 
(£,i(t)^j(t')) = 2TSij6[t — t']. Working out this expression gives 



±6 



cos(A ik ) sm(6i-6 k ) + sin(A ik ) cos(0j-0 fe ) 



-hemffli- (3) 

with Aik = Aik for i < j, and A^ = —Aki for i > j. We introduce disorder by 
drawing the pinning angles <pi independently at random from a distribution p{<j>), 
and the relative angles A^ independently at random from a distribution P(A), with 
J dA P{A) cos(A) = (to ensure that the bond disorder will retain significance 
for N — > oo). This system will generally exhibit competition between alignment 
of the spins to the pinning fields and the realization of prescribed relative angles 
between pairs, and have a high degree of frustration. Each pair of spins try 

to achieve 0j — 6j + A\j = 0, so a bond-frustrated loop {ii,...,^} is one where 
Ai t i 2 + Ai 2 i 3 + . . . + A iL _ r i L + Ai Lil 7^ (mod 2tt). We simplify the bond average 
by choosing the {Aij} to be binary: P(A) = ±5[A-A*] + \8[A-A*--k). Without loss 
of generality we may take — 5 < A* < ~. However, only for A* £ {0, 4-7!"} will it be 
possible to calculate the disorder-averaged free energy for the system (H) in terms of 
a standard replica mean field theory; these two cases are the subject of our paper. 
Our first (conventional) model corresponds to P(Aij) = ^S[Aij] + i5[Ay — n]: 

Model I: H = — -= ^ J i3 - cos(^-^) - /i^cos(0<-<&) (4) 

i<j i 

with Jij = cos(Aij) G { — 1,1}. In this model energy minimization will translate for 
any pair (i,j) of oscillators into the objectives 

Jij = l: 0i-Oj^O (mod 2tt) 

Jij = — 1 : di — 9j^ir (mod 2tt) 

(i.e. ferro- and anti-ferromagnetic interactions, respectively). The Hamiltonian (^) 
defines a Gibbs measure corresponding to the Langevin forces /, = —dH/dOi = 
— (2K/^/N) J2k^i sin(9i—0k) — hsin(8i—(f)i), provided we define = Jjj. A triplet 
of spins (i, j, k) is now bond-frustrated if an odd number of the relative angles involved 
equals —1, and unfrustrated otherwise (see figure [l]). Our second (less conventional) 
model is obtained for P(Aij) = ^[Ay-f] + \5[A i:j + ^]: 

2K 

Model II: H = — -= ^ J i4 sin^-fl.,) - ft^ cos(0;-<^) (5) 

i<j i 
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Figure 1. The low energy states of the two models (kj) (upper row) and (n) 
(lower row), for TV £ {2,3} and in the absence of pinning fields. Arrows indicate 
the orientations of the vectors Si = (cos((9i), sin(9i)). Without pinning fields only 
relative angles are important, so we may put 8\ = ^n. The bonds Jij £ {—1, 1} 
are drawn as solid (1) and dotted (—1) line segments, respectively, which connect 
the vectors Si and Sj. Frustrated clusters, where no state exists which minimizes 
all terms in the Hamiltonian simultaneously, arc indicated with '?'. In model I the 
triplet (1,2,3) is bond-frustrated if an odd number of the bonds { J12, J13, J23} 
equals —1. In model II the triplet (1, 2, 3) is always bond-frustrated. 



with Jij = — sin(Ay ■) S { — 1, 1}. Here energy minimization will translate for any pair 
of oscillators into the objectives 

Jij = 1 : Oi-63 -> tt/2 (mod 2tt) 

Jij = -1 : Bi-6j -> -tt/2 (mod 2?r) 

The Hamiltonian (^) defines a Gibbs measure corresponding to the Langevin forces 
/, = —dH/dOi = —(2K/y/N) J2k^i Jij cos(di—9k) — hsin(di — <f>i), provided we define 
Jij = —Jji. Now every triplet (i,j,k) is bond-frustrated, since with Ay G {— 5, 5 } 
one always has Aij + Ajk + Aik = £n/2, with £ integer and odd (see figure 0). 

It is not a priori clear whether and how the solutions of models (Q) and (||) are 
related. For instance, finding their ground states for N — 3 reduces to minimizing 
Hi{9 1 ,d 2 ) = Jicos(0i) + J 2 cos(# 2 ) + J 3 cos(d 1 -6 2 ) and H u (e 1 ,6 2 ) = Jisin(^) + 
J2sin(#2) + J3 sin(^i — 62), respectively (since rotation invariance allows us to put 
#3 = 0), with Ji E { — 1,1} (chosen randomly). Solving this simple problem reveals 
that the ground state energy E\ of H\ is dependent on the bonds { Ji}, E\ S {—1,-3}. 
Averaging over the bond realizations gives E\ = —2. In contrast, for the ground state 
energy of model II one finds En = —1^3, independent of the bond realization. 

Pinning fields break the global rotation invariance of (||) which is present for h = 0, 
and simplifications of our equations will thus have to be based on symmetries of the 
pinning field distribution. Any non-zero value of the average pinning direction (cj)}^ 
can be transformed away, via 0i — * + {4>)$, without affecting the interaction terms 
in our Hamiltonian. Thus, without loss of generality we may choose (<p)cj, — 0. In the 
remainder of this paper we restrict ourselves, for simplicity, to reflection-symmetric 
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pinning field distributions, where p{<f>) — p(~4>) for all G [— tt, tt]. However, the free 
energies of our models (^,^) are both invariant^ under the transformation 



i 



(l+7])p(<f) - A 



7T . 1 , 



(6) 



for any 77 e [—1,1] and A € [0, 27r], so we solve implicitly for a larger family of 
models. For high temperatures we expect the symmetries of p{4>) to be inherited 
by the solutions of our models. We define a complex (single-site) disorder-averaged 
susceptibility x, upon adding small perturbations to the external pinning angles (j)^. 



X 



-y 



iNh 



-i(pi 



d 



(7) 



Here 



and [. . .] denote equilibrium and disorder averages, respectively. For 



perfectly linear single-site response, i.e. (e 1 ' 



x e 



definition (|7]) would give 



x/h. Evaluating (Q) for systems with Hamiltonians of the form (||) gives 



X = 



iN ^ 



sin 



(8) 



3. Equilibrium analysis 

3.1. Solution via replica theory 



We calculate the disorder-averaged free energy per spin / — — lim7v^ 00 (/3A r ) _1 log Z 
for the Hamiltonian (|^), using the identity log Z = lim n _,o log Z n and the standard 
manipulations of conventional replica theory p9| , giving 



/ = — lim lim 



1 



AT->oo n->0 f3Nn 



log 



U de ° 



.a = l 



The disorder average, with P{A) = ±5[A-A*] + \S[A- A* -it], gives: 



i<j 



N 



-A*)} 



Expansion for large JV, however, only leads to an expression in terms of the usual 
single-site replica order parameters if sin(2A*) = 0, i.e. for A* e {0, §} (f||||). Here 



with 



i i 



t. The transformation (0) implies <j>i —►</>.; + A + t; ? , where Tj E { — 1, !}• The variables A and {r;} 



can be gauged away by putting 



+ \ + Ti^ in combination with Jij — > Jjj cos(J (Tj —75)). 
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We proceed as usual, isolating observables via appropriate ^-distributions (which 
introduces conjugate order parameters), leading to site factorisation. The result is 

7 = - lim \ cxtr {$[{§**, q**}] + *[{«**}]} (9) 

i 



*[■■■] = a E E + ^ 2 cos2 (^) E E^] 2 

** a/3 ** a/3 

+ 2/3X 2 sin 2 (A*) ]T [C/3C/3 - C/3C/3] (10) 

a/3 

/?*[...] = (log y"de M(e|{9**})>* (11) 

with (**) e {(cc), (ss), (sc), (cs)}, and with 

M(0|{q**}) = e f3h ^-"> cos ^ e °^- l_i Sc,fi^S cos ( e a) cos ( e /3)+9a0 cos(e Q )sin(e /3 )] 

x e - i I] c ,^K s /3 sin ( 9 °) sin (^)+«S sin ( e °) cos (^)] ( 12 ) 

One always has q^ s = q s p c x , which allows us to simplify the saddle point equations. For 
both models we find upon varying the conjugate order parameters {q^}'- 

_ / Jd9 coB(0 x )cos(9 p )M(O\{q**}) \ 
Qx " \ JdO M(0\{q**}) J, [i6) 

ss / Jd6 sin(fl A )sin(g p )M(0|{g**}) \ 
9a » \ fMM(6\{q**}) U ( > 

sc / fdB S m(0 A )co8(e p )M(g|{q**}) \ 

^ \ fd9M(e\{?*}) I, 

The models I (A* = 0) and II (A* = f ) differ only in the equations resulting from 
variation of the order parameters {9*^}: 

Model I: (16) 

Model I • ^ = = ( i 71 

' = -2i0*lCq%, = -2i(3 2 K 2 q% { } 

3.2. Replica- symmetric solutions 

We now make the replica-symmetric (or ergodic) ansatz (RS) for the saddle-points: 

Qa*0 = (Q**-q**)S a p + q±± (18) 

q% = U(2pKf(R^-r^)5 afj + U{2pK) 2 r^ (19) 

Insertion into the saddle-point equations tells us that 

R cc = cos 2 {A*)Q cc + sin 2 {A*)Q SS r cc = cos 2 (A*)q cc + sin 2 {A*)q ss 
R ss = cos 2 (A*)Q ss +sin 2 (A*)Q cc r ss = cos 2 (A*)q ss +sm 2 (A*)q cc 
R cs = cos(2A*)Q cs r cs = cos(2A*)q cs 
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The identity Q cc + Q ss — 1 allows us to reduce the problem to solving five coupled 
non-linear equations. We will compactify our notation by using the short-hands 

JdO f(9)M(6\x,y,u,v,<t>) 



= JdxDvDuDv... (f(0)U = 



Jd9 M{6\x,y,u,v,4>) 
We can write the five remaining saddle-point equations as 

Q cc = «( (cos 2 (0)), ))) Q cs = ((( (cos(0) Bin(0)>* )))* (20) 
q cc = ((( (cos(e))l q ss = ((( (sin(0)) 2 »), (21) 

q cs = ((( (cos(0))*<sm(0))* »>* (22) 
Tracing back the physical meaning of these observables gives: 

Qcc= hm -i V (cos 2 (00), Qc,= Um i- V (cos(^) sin(^)> (23) 

iV— >-oo iv z — ' N—*oo iV ^ — ' 

9cc = Hm l^cos^)) 2 , Qss = hm l^< s in(^)) 2 (24) 
9cs = hm 4^< cos (^))(sin(^)) (25) 

i 

Our two models differ only in the form taken by the effective single spin measure 
M(9\x, y, u, v, <p) in the RS saddle-point equations: 

Mi(6\x, y, U, V, <j>) = e P hc °< ( H>) + (l 3K ) 2 [[ 2 Q^- 1 +'!^-'!^ co S (2e)+2(Q cs - 9cs ) sin(20)] 
J3K cos 

Ml(6»|x, U, U, (f>) = e / 3 ' l ™ s ('M)-(/3i<") 2 [[2Qec-l+9 3s -gcc] cos(20)+2(Q cs - gcs )sm(29)] 

x e /3Kcos(e)V2[xV2^^+(iu+v)^^+pK S in(e)V2[yV2^^+(iu-v)^q^] ^7) 

(apart, in both cases, from a constant pre-factor e^ K "> ^-i^-q^Ay Xo compactify 
future notation we define the following short-hands 

7cc = (cos 2 (0)K-(cos(#)) 2 lss = <sin 2 (fl))*-(sin(0)) 2 (28) 

7cs = (cos(0) sin(0))*-<cos(0))*(sin(0))* (29) 

Note that in replica-symmetric states one has (((7cc)))</> = Qcc — Ice, (((7ss)))^ = 
1-Q CC -<? SS , «(7cs)))^ = Qcs-Qcs, and |7 CS | < ±(7 CC + 7 SS ). 

Similarly we can work out the RS free energy per oscillator, by insertion of the 
RS ansatz into (^), followed by standard manipulations. This results in 

/ = extr^^/^Q**,^*}] (30) 
/[. . .] = (3K 2 U - ^ J DxDyDuDv log J d9 M(9\x, y, u, v, (p)J (31) 

with the effective measures M(9\ . . .) defined in (||J27|), and with 

Ui = 2Q cc (Q cc -l) + 2Q 2 CS - q cc {q cc -l) - q ss {q ss -l) - 2q 2 cs (32) 
Un = - 2Q cc (Q cc ~l) - 1 - 2Q 2 CS + q cc + q ss - 2q cc q ss + 2q 2 cs (33) 
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The nature of the extremum in ( p0[ ) follows from the high-temperature state. 
Expanding for (3 -> gives /[. . .] =-i"log(27r)-i/?/i 2 + 2f3K 2 exii{tj} + 0{(3 2 ) 1 with 

1 



Ul =Q 2 cc-Qcc + Ql s - T:(q 2 1 



ice ' Qss 



Qlc + Qc 



l 



q cc q s 



For (3 — > we must find the paramagnetic solution Q cc = ^ and Q cs = q cc = q ss = 
q cs = 0, so the nature of the desired extremum in (J3TJ) is 

I: min w.r.t. {Qcc, Qcs}, max w.r.t. {q cc , q SSl q cs } (34) 
I: min w.r.t. {q cs , q C c~ q ss }, max w.r.t. {Q cc , Q cs , 

Qcc ~r Qss } (35) 

Finally, for RS solutions we can also work out expression (|^) for \ in the limit N — > oo: 

^ = «( (sm 2 (0-0)K-<sin(tf-0)) 2 »>, 

-i «( (Bin(fl-0)oos(d-^))*-(Bin(d-0))*(cos(fl-^))* »}* 
Equivalently, with the conventions (^8|,^9]) this can be written as 

= 7}{l-qcc-qss) - ((( icos(20)(7 cc -7 ss ) +sin(20)7 cs 

- t ((( cos(20)7 cs - -sin(2(^)(7 cc -7 ss ) (36) 



3.3. The AT lines 

Here we determine the stability of our RS solution against 'replicon' perturbations: 

1a*/3 9^ RS + T Q with C/? = Vfa and C/j = ( f ° r &U ^« = for &U 

X] a J)*]j = Sal^a = for all a, and with all |?7*^| <C 1. Details of the derivation 
of the equations signaling the this type of instability (the AT [^8| lines), are given 
in Appendix B. Our calculation is more complicated than that in e.g. 123, E7J, 



because, due to the pinning fields, we cannot generally use rotational symmetry; 
replicon fluctuations can now have a more complicated structure. We found two 
types of AT instabilities (the physical RSB transition associated with these is the one 
occurring at the highest temperature). The first is the same for our two models: 

Models I & I : (T/2A) 2 = (((icdss -ll)))^ (37) 
The second AT instability is found to be model dependent: 



with 



Model I : 
Model I 



C = 



Det[E - (T/2A) 2 !] = 
DetpS - (T/2A) 2 C] = 




mii 2 cc+' 



\ 2 V Ice is 

(((7cs(7cc+7 



7c 2 s »>0 {i\{l 2 cc-l 2 s))))0 («7c S (7ee + 7 SS )}»0 

^ (((|(7 c 2 c + 7 s 2 s )-7c 2 s »>0 (((Jcshcc-lss)))), 

}<f (((lcs{7cc-Jss))))<p (((lccjss+lcs)))<l> 



(38) 
(39) 



(40) 



(41) 
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4. States with global reflection symmetry 

We inspect the effect of global reflection — > —di, on the order parameters (within 
the RS ansatz), with the help of the identifications ( p3| , p4p5| ) : 

9'i = -0 t for all i: ^ Z 1^"' , _ ( 42 ) 

Ice — 9cci Qcs — ^csj Qss — Qss 

Invariance under this transformation implies Q cs = q cs = 0, and breaking of the 
associated symmetry is signaled by a bifurcation of Q cs ^ and/or q cs 0. 

4-1- Implications for free energy and order parameters 

For reflection-symmetric states the Gaussian variables {u, v} disappear from the 
problem, and the effective single-spin measures ( p6| , p7| ) simplify to 

Mi(9\x, y, 4>) = e / 5/lcos (e-0)+(/3J<") 2 [2Qcc-i+g ss -gcc] cos(2e)+2pK[ C os(e)x^q^+sm(6)v^g7^ 

(43) 

M^(6\x,y, (f) = e l 3hcos (^4>)-(l3K) 2 [2Q cc -l+q ss -q cc ]cos(29)+2[3K[co S (9)x^^+sin(e)y^] 

(44) 

These expressions show, in combination with the underlying symmetry p(cj>) = p(—(f>), 
that for any set of functions {kg} (with I = 0, 1, 2, . . .): 

((( k {-4>) J] vM-0)>* »>0 = ((( *o(0 \{{Ue)), ))h (45) 

i>0 £>0 

One is left with just three coupled order parameter equations: 

Q cc = ««cos 2 (0)K)» q cc = ««co S (0))*»V q ss = ««sm(0))*»U(46) 
Inserting Q cs = g cs = into ( [SOp^ , |33| ) shows that for reflection-symmetric states the 
free energy per oscillator equals / = extr { Q cc (Zss , qcc }f[Q cc , q ss , Qcc], with 

fi[Qcc,q ss ,qcc\ = PK 2 [qcc+qss+ZQcdQcc-^-qic-ql 

1 

fa[Qcc,q S s,qcc] = PK 2 [q cc + q ss -2Q cc (Q cc -l)-2q ss q cc -l) 



DxDy log / d6 Mi(6\x,y, 



^jDxDy log J d6 M R (e\x,y, 



From these expressions one can easily extract the replica symmetric ground states, 
using a saddle-point argument as [3 — ► oo in the entropic term: 

lim fi[. . .}//3K 2 = q cc +q ss + 2Q cc (Q cc -l)- q 2 c - q 2 ss -\2Q cc -l+q ss - q cc \ 

p^oo 

lim /n[. . ]/(3K 2 = q cc +q ss -2Q cc (Q cc -l)-2q ss q cc ~l- \2Q cc -l + q ss -q cc \ 

According to (|34||55|), we must for model I minimize with respect to Q cc and maximize 
with respect to {q cc ,q ss \, whereas for model II we must minimize with respect to 
qcc~~ qss and maximize with respect to {Q cc , qcc+qss}- For both models this leads to 

qcc+q ss = l, Qcc = q cc - q ss ) (47) 

giving in both cases extrQ cci9ssj(?cc lim /3 ^ 00 (/3i ; ir 2 ) _1 /[(3 cc^ss^qcc) — 0- Hence there 
is a finite ground state energy. To find the RS ground states of our two models from 
the family q cc -\-q S s = 1, one would have to inspect the next order in T. In addition 
we must, of course, expect replica symmetry to be broken for low temperatures. 
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Reflection symmetry is also found to simplify expressions (|3q , |39| ) for the AT 
instabilities, due to (((7cc7cs))) <t> — (((TssTcs)))^ = 0. In particular, with (filf ) one 
can simplify the matrix (|4l]) to 

(((^(7 C 2 C +7 S 2 J+7 C 2 S )»0 («i(7c 2 c-7 2 s )>»* 
{(\illc-lls))))* (((|(7 C 2 C +7 2 S )-7 C 2 S )»0 

(((7 CC %,+7 C 2 S »>0 

In combination with (|37|), and upon rejecting solutions which are immediately seen 
not to give the highest transition temperature, the AT lines of our two models are 
found to be the solutions of the following equations, respectively: 

2 



(48) 



I 



I 



T 
2K 



T 



1 



= max 



1 



<«7 2 C +7L)))0 + V(((7 C 2 S )»^ + 7(((7 2 C -7 2 S )))^ , (((7 C c7-+7 2 s »>* 



2K 



max {(((7 2 s »>0 + V / (((7c 2 c )»0(((7 2 s »)0 , (((icdss 



(49) 
(50) 



For both models the left of the two arguments in the corresponding extremization 
problems give the required maximum. To see this, one first notes that 



«(7cc7 SS »>0 = \{{{llc +7 S 2 S - (7cc-7ss) 2 ))U < \(((l 2 cc+l 2 ss)))* 



(51) 



1 



1 



<«7 2 c )»0<«7 2 ,»)0 = g {<«7 C 2 C +7 2 S »>^ ~ «(7 C 2 C »)^ ~ <«7 2 c»> 2 } > ^H c +ll)))l (52) 
For model I we subtract the right argument of (J49) from the left argument and find 



LA — RA = I(«( 7cc - 7ss )2))) - <« 7 f s ))), + ^{{{ 7 l)))l + \{{{ 7 l- 7 l)))l 



> l((((lcc-lss) 2 )))^ > 



Thus the maximum in (^) is always realized by the left argument. Similarly, upon 
subtracting the right argument from the left argument in (^0|), using (51 152]), we find 

LA - RA = 2((( 7c 2 J)), + y / itM^- («7«7«>»* 

> ^{Uc)))l{hl)))l- = ^V2-l)((( 7c 2 c + 7 L>»0 > o 



Hence also the maximum in ([30|) is realized by the left argument. We conclude that 

(53) 



AT line I : 
AT line E 



T 



2K 

T 
2K 



J(((7 2 c +7L>>>0 + V(((7c 2 s >>> 2 , + l(((7 2 c -7 2 s >)> 2 , 



= (((7 2 S »)0 + V(((7 2 C »>0(((7 2 S »>0 



(54) 



Finally, in reflection symmetric states with p(4>) = p{—<j>) V0, where we may use 
identity (fi"5"|), the replica-symmetric susceptibility ( p6|) is found to be purely real: 

Xrs 1 1 

— = ^-q cc -q ss ) - (((- cos(20)( 7cc - 7ss )+sin(20) 7cs ))) (55) 
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4-3. Reflection symmetry breaking transitions 



Assuming reflection symmetry breaking (to states with Q cs =/= and/or q cs ^ 0) to 
happen via second order transitions, allows us to determine its occurrence by studying 
the relevant entries of the RS Hessian of the free energy. Such transitions occur when 



Dct 



d 2 f[Q cc , q33 , q J/dQl 

d 2 f{Qcc,q ss ,q cc }/dQ CS 6 



d 2 f[Qcc,q BB ,q cc }/dQ cs C 
d 2 f[Qcc,q ss ,qJ/dq 2 c 



(56) 



The relevant matrix elements of the Hessian are given in Appendix A . We can deal 
with both our two models simultaneously upon defining the variable r G { — 1,1}, 
where r = 1 for model I and r = — 1 for model II. This allows us to write: 

d 2 f 



1 

m 2 dQ 2 cs 

1 d 2 f 

PK* dq 2 cs 
1 d 2 f 

f3K 2 dQ cs dq c 



At - 8((3Ky\ 1 



4t + 8{(3K) 2 \ 2 



They involve 



Ax =-<<((sin 2 (20)K-<sin(20))2))) 



(57) 
(58) 
(59) 

(60) 



A 2 = ((([(sin 2 (0))^(sin(0)) 2 ][(cos 2 (0)),-3(cos(0)) 2 ] 

+ [(cos 2 (0)),-(cos(0)) 2 ][(sin 2 (0)),-3(sin(0)) 2 ] 
+ 2[(sin(0) cos(0))*-(sin(0))*(cos(0))*][(sin(0) cos(0))*-3(sin(0))*(cos(0))*])))^ 

(61) 

A 3 = 2<(«sin(20) cos(0))*(sin(0))* + (sin(20) sin(0))*<cos(0))* 

- 2<sin(20))*<inn(0)}*<a»(0)}*}}>* (62) 
Insertion of the second derivatives into ( |56| ) , using the above short-hands ( |60| , |6l| , |62] ) 
reveals that reflection symmetry breaking transitions are marked by the highest 
temperature for which the following functions Ejjj(T) become negative: 



Model I : 
Model II 



Zf(T) = (T/K) 2 — V /(A 1 -A 2 )2-A§-A 1 - A 2 



Eff f (T) = (T/K) 2 - J(A 1 -A 2 )2-A| + A : + A 2 



(63) 
(64) 



Since the {A^} are bounded, expressions ( |63]]64| ) confirm that reflection symmetry will 
always be stable for sufficiently high temperatures. 



5. States with global rotation symmetry 

For uniformly distributed pinning angles, i.e. p(4>) = (27r) -1 , we may expect the 
macroscopic state to have global rotation symmetry for sufficiently high temperatures, 
in addition to reflection symmetry. We therefore inspect the effect of global rotations 
0i — > Oi+ip on the order parameters (within the RS ansatz), with the help of (23 2^ , |25| ): 

8'i = &i + 4> forall i: q' cc + q' ss = q cc + q ss (65) 
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Q'cc-l\ = /cos(2^) -sin(2V)\ (Qco-\ 
Q' cs J { sin(2^) cos(2^) J \ Q cs 

Wcc-Qss)\ = fcos(2V>) -sin(2V)\ ( \{q cc -q s 
q' cs ) \ sin(2^) cos(2^) J \ q cs 

Invariance under all global rotations implies Q cc = i, Q cs = q cs = 0, q cc = q ss = q 
(rotation- invariance implies reflection- invariance), leaving just one order parameter in 
rotation invariant states. The invariant manifolds in order parameter space are 

{Qcc~\? + Q 2 cs =ei (66) 

^(q cc -q ss ) 2 + q 2 cs = e 2 (67) 

with the invariant state corresponding to ei = £2 = 0. We will also find instances 
of rotation-invariant states without the pinning field distribution having rotational 
symmetry. Insertion of {Q cc — h, Qcs — qcs — 0, q cc = q ss = q] as an ansatz into our 
RS order parameter equations reveals that (for nonzero h) the following condition is 
necessary and sufficient for the existence of a rotation-invariant solution: 

(cos 2 (0)) = (sin 2 (0)) (68) 

Hence, unless explicitly stated we will in this section not assume p{4>) to be uniform, 
but rely only on the two properties p(tfi) = p(—<fi) (assumed to hold throughout this 
paper) and (cos(2</>))0 = (to guarantee the existence of a rotation-invariant state). 

5.1. Implications for free energy and order parameters 

For rotation-symmetric states the measures ( ^3||4i| ) become identical: 

M(9\x,y 7 4>) — e hcoB ( e -<l>)+ 2 P K Vql xcos ( e )+y sin ( e )] (69) 

For any value of <f> one can carry out a suitable rotation of the Gaussian variables 
(x,y) to eliminate </> from the measure (|69[), leading to the following general identity 
for any set of functions {hi} (with I = 0, 1, 2, . . .): 

e>o £>o 
where (. . .) refers to averages calculated with the ^-independent measure M(0\x, y) — 
M{6\x, y, 0). As a consequence our subsequent calculations for rotation-invariant 
states will repeatedly involve various derivatives of the following generating function: 

Z[x,y] = log J dO e /^ c ° s ( e )+ 2 /3^V9[zcos(e)+ysin(0)] 

= log(27r) + log/ [S] (71) 

with the short-hand 5 = f3y/ (h + 2Kx^/q) 2 + (flKy^fq) 2 , and in which I n [z] denotes 
the n-th modified Bessel function |30]. For instance: 

Wkjlfcl*>v] =<sm(0)> o (73) 

' <>2 Z[x,y] =(cos 2 (0)) o -(cos(0))2 (74) 



{2(3K^qf dx 2 
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1 d 2 



(2(3K^q) 2 8y 2 

i d 2 



Z[x,y] = (sin 2 (#)) -<sin(#)) 2 (75) 
Z[x,y] = (sin(0)cos(0)) o -(sin(0))o(cos(0))o (76) 



(2/3K y /q) 2 dxdy 
1 d 3 

(2f3K V qrdx 2 dy Z[x ' y] = (^W cos 2 (0)) o - (sin(0)) o <cos 2 (0)) 
1 d 3 



(2(3K y /q) 3 dxdy 2 



-(sin(2#))o(cos(0))o + 2(sin(0))o(cos(0)) 2 (77) 
;Z[x,y] = (sin 2 (0) cos(6»)) -(cos(6»)) (sin 2 (6i)) 

- (sin(20)) o (sm(0)) o + 2(cos(0)) o (sin(#)) 2 (78) 

In order to work out the remaining order parameter equation for q we apply the 
identity @ to ki(6) = k 2 (0) = cos(0) and to fei(0) = k 2 (0) = sin(0), giving 

((((cos(0)) 2 ))) = (cos 2 (0)) (( (cos(0)) 2 » + ( S in 2 (0)) ((( S in(0)) 2 )) (79) 
({{{sin(d))l))U = (sk 2 (0)>4 (cos(0)) 2 )) + (cos 2 (0))4(sin(0)) 2 » (80) 
Thus the remaining order parameter q is, for both models, to be solved from 
q = l - J DxDy ( (cos(tf)) 2 + (sin(tf)) 2 ) 

= i y DxDy [(cos(0)) 2 + <sin(0)) 2 ] (81) 

Using properties of modified Bessel functions, such as -^I n [z\ = \ (In+i [ z ] + -^n-i[-z]) 
and I n [z] — (z/2n)(I n -i[z] — I n +i[z]), the two relevant expressions ( |72] , |73| ) give 

^.- ""^ 'M (82) 

(sin(0)) o = — py (83) 

After writing (x, y) in polar coordinates, one then finds ( pl| ) reducing to 

with 



3(r, VO = (3^4:qK 2 r 2 + h 2 + 4K^r^cos(^) (85) 

Inserting Q cs = q cs = 0, Q cc = |, and q cc — q ss = q into expressions (|3(],|3^,|3|) 
shows that for rotation-symmetric states and for both models I and II the disorder- 
averaged free energy per oscillator equals / = extr g /[g] — 4 log(27r), with 

f[q] = -2(3K 2 {q-\) 2 - - [**t f°dr re'^ logJ [H(r, V)] (86) 
2 P Jo I" Jo 

According to ( [34|]35| ) the function f[q] is in both cases to be maximized, so the replica 
symmetric rotation-invariant ground state would be q = i. 
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Various terms involving the covariances 7** (£8|,g9|) can be simplified via (|70|). We 
define 7^ = 7**10=0 (i.e. as calculated with the <p = averages (. . .) ). For instance, 
focusing on the terms occurring in and using (|68|): 



(«7 C 2 C >»0 = <(|(7 CC +7 S J 2 + 7L» + <cos 2 (20))0((i(7 cc -7 ss ) 2 -7L 

(«7L)))0 = <(J(7 CC +7 S J 2 + 7L» + <cos 2 (20))0((i(7 cc -7, s ) 2 -7 2 s 
,1 



(87) 



= <sin 2 (20)) «( i (7 cc -7 s J 2 )) + <cos 2 (20)) «7 2 s )) 

|«((7cc-7 s ,) 2 »)0 = <sin 2 (20)) ((7D) + <cos 2 (20)) <(i(7 cc -7, s ) 2 » 
From ( |7"4] , f75| , f76| ) , in turn, one obtains explicit expressions for the required {7**}: 



7 CC 2 

_ 1 

7 SS - 2 



l3 2 {h + 2K Xy /qY 



m 

2I [E] ' E 2 
2/o [S] + 



m 

Jo [HI 



^2 



7 C 



(3 2 {h+2Kx s /q){2Ky y fq) 



—2 



Jo [3 



(89) 
(90) 

(91) 
(92) 
(93) 



Combination of the above results with (|53|,|3J) gives the equations for the AT lines. We 
note that for rotation invariant states the two expressions (^,Q) become identical: 
1 1 2 



4((^7 2 C + ^7L+7 2 S 



/2IHI 



= / DxDy ■ 



1 - 



After transformation of the Gaussian variables to polar coordinates, this gives for the 
AT lines of our two models the appealing result 



~ T~ 


2 




- 


K 









/ a [S(r,V)] A 2 [S(r,V0] 



/o[S(r,V)] / 2 [S(r,^)] 



1 



(94) 



In a similar manner we work out expression ( pq ) for rotation-invariant states: 
^ = -(l-2g) - <«-cos(20)(7 cc -7 ss ) + sin(2,/))7 cs )))0 



i(i-2 9 )-i (7, 



2^ 



ar re 2 



7 SS ) 

/ 2 |H(r,V)] A 2 [2(r,^)] 



7 [S(r,V)] I 2 [S(r,V)] 



r h 2 +4K 2 r 2 q cos(2ip)+4hKr,/qcos(ip) 



(95) 



/i 2 + 4if 2 r 2 g + 4/i-KV cos(^) 
with the abbreviation (pq). Note that lim^o S(r, 1/') = ^PKr^Jq (i.e. independent of 
?/;), so that for rotation invariant states one has 

lim xrs = P{\ - q) (96) 

ft— »0 z 
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5.3. Rotation symmetry breaking transitions 

According to ( |66| , |67| ) we can study rotation symmetry-breaking bifurcations most 
conveniently after the following transformation of the order parameters: 

Qcc = g + £ i cos(wi), Q cs — eisin(wi) 

q cc = q + £2 cos(w 2 ), q ss = q - e 2 cos(w 2 ), q cs = £2 sin(w 2 ) 



Reflection symmetry breaking transitions are already marked by ( p3|]64| ) , so we 
restrict ourselves to rotation symmetry breaking transitions which preserve reflection 
symmetry, i.e. Q cs = q cs = 0. Thus we put u>\ = u>2 = and focus on 

Q cc =2 + e i> Qc S = q cs = 0, q cc = q + £2, q ss =q-^2 

Our interest is in continuous bifurcations of e\ ^ and/or e 2 ^ 0. We expand the 
free energy per oscillator /[...] to be extremized (|30|) around t he rotation invariant 
solution, using the second order derivatives in Appendix A. Due to reflection 
symmetry, only those second derivatives with an even total number of 's' subscripts in 
the two corresponding order parameters can be nonzero. This gives for <C e±, e 2 <C 1 
(with /ri[. . .] denoting the free energy of the rotation invariant state): 

/[..., _ /Rl[ ..., i ? /|V + !V 2 av 



2 1 dQ 2 cc " 2 2 \ dq 2 cc 8ql 8q cc dq ss 

I OQccOqcc OQccOqss J 

Hence the continuous rotation symmetry breaking transitions are marked by 

d 2 f/dQ 2 cc d 2 f/dQ cc dq cc -d 2 f/dQ cc dq ss 

d 2 f/dQ cc dq cc - d 2 f/dQ cc dq ss d 2 f/dq 2 cc + d 2 f/dq 2 s - 2d 2 f / dq cc dq S! 



Det 



(97) 



We work out the second derivatives for rotation invariant states (where the effective 
measures of the two models are identical), using (f70|). Given our requirements 
p(4>) = p{—4>) and (cos(20))^, these are found to depend on the choice made for the 
pinning field distribution only through (cos{A<f>)) Again we deal with both models 
simultaneously upon defining r G { — 1,1}, where r = 1 for model I and r — — 1 for 
model II. This allows us to write: 
1 d 2 f 



= 4r - 8((3K) 2 Pl (98) 



(3K 2 dQ 2 cc 
1 (d 2 f d 2 f n d 2 f 



(k-w* 2 ps7r- 4T+W w (99) 

where straightforward but tedious bookkeeping shows the quantities {pi, p 2 , p^} to be 
Pi = - (cos(20) 2 - (sin(20) 2 )) + i(cos(40)) (((cos(40)) o - (cos(20) 2 + (sin(2fl) 2 )) 

(101) 

P2 = (( [(cos 2 (0)) o - (cos(0)) 2 ][(cos 2 (#)) o - 3(cos(0)) 2 ] 
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+ [(sin 2 (tf)) - (sin(0)) 2 ][(sin 2 (tf)) o - 3<sin(0)>2] 
-i[<cos(20)) o - (cos(#)) 2 + (sm(0))l][(cos(29)) o -3(co S (e))l+3(sm(9)) 2 o ] )) 
1 , 
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(cos(#))4 ^[(cos(20)) o -(cos(0)) 2 



sin(0))2][(cos(20)) o - 3<cos(0))2+ 3<sin(0))2] 

-2[(sin(0) cos(0)) o - (sin(0)) o (cos(0))„] [<sin(0) cos(0)) o - 3(sin(6»)) (cos(6»)) ] )) 

(102) 

p 3 = 2(((sin 2 (0)) o (cos(0)) 2 + (cos 2 (9)) (sm(9)) 2 - (sin(20)) o (sin(0)) o <cos(0)) o )) 

+2(cos(4^)) « (sin 2 (#)) (cos(fl)) 2 + (cos 2 (#)) o (sin(0)) 2 + (sin(20)) o (sin(0)) o (cos(0)) c 
- 2(sin 2 (6»)cos(6»))o(cos(e')) - 2(cos 2 (6>) sin(0)) o (sin(6»)> » (103) 



Insertion of the second derivatives into (p7\), using ( |10l]Jl02|J103| ), shows that rotation 
symmetry breaking transitions which preserve reflection symmetry are marked by the 
highest temperature for which the following functions become negative: 



Model I : 
Model II 



Sj (T) = (T/K) 2 - J( P1 -P2) 2 -P 2 3 -PX-P2 



(104) 
(105) 



EJf(T) = (T/K) 2 - ^( Pl - P 2) 2 -Pl + Pi+P2 

What remains is to work out (101,102,103). For most terms we can use (71 7^-F^,^- 
P3|). One further object in pi, (((cos(40)) o )), can be calculated directly as follows 



«<cos(40)) o » = (( 



JdO cos(40 + 4atan[ 7r f§^])e 



27rJ (S) 
; ( cos(4 ^n\ 2Ky T f ])^g 



1 - 



h + 2Kx^/q' ' I [3]" 
4 (2^yV9) 2 (^ + 2 ^V9) 2 



The final result of the exercise is, after transformation of (x, y) to polar coordinates: 



dr re ? r / — < 1 

7T 



(cos(4(/))) 



o tt Uo[S(r,V)] / 2 [S(r,V)] 
8/3 4 [2Xr^g sin(V>)] 2 2Kr^/q cos(V>)] 2 

1 s^MO 



(106) 



P2 



2 Jo 7o TT 1 I 2 [S(r,^)]i I / 2 [S(r, 



(cos (40))^ / dr re 2 



i r2 r# f/ 2 [3(r,V)] iflHM)]] f/ 2 [HM)] 3/ 2 [S(r,V0] 



o tt Uo[S(r,V)] / 2 [S(^)]j(l„[H(^)] / 2 [S(r,^)] 
8/3 4 [2^rVgsin(^)] 2 [fe+2grygco S (^)] 2 | 
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7o[H(r,V0]J ItU^)] 

ff #f I 2 [Z(r,i>)}\ ( If[E(r,4>)] 2J 2 [H(r,^)] 



Jo|H(r s V)]J Uo 2 [S(r,V)] /o[S(r,^)] 
8/3 4 pKr^/g sin(^)] 2 2*0^/3 cos(V')] 



(108) 



The general identity J n+ i[z] = J„_i[z] + nl n [z] (I\ [z] — I [z] ) jl\ [z] allows us to 
express any I n >2[z] in terms of the trio {ioN) h[z], h[z]}- Since the {p^} are 
bounded, expressions (104, |105| ) confirm that, if rotation invariant states exist (i.e. if 
(cos(2(^))0 = 0), rotation symmetry will be stable for sufficiently large temperatures. 

To round off our discussion we return to reflection symmetry breaking transitions. 
We work out (^,^l|,^2|) for rotation invariant states, and find (after a lengthy but 
straightforward calculation) that 



Al = 4 



do Jo - V /„ 2 [S(r»]/ 



4^"™^ - Jo tt Uo[S(r,^)] I 2 [S(r»] 

SP 4 [2Kr y /qsin(^)} 2 [h+2Kr v /qcos(^)} 
S 4 (r,V) 



_^ u/-j i yij iHT^uy^uj^j] ^ (109) 



■kcntm* Hair re~^ f ^ ( Mg^ - ZL^^ l ( Mg^ _ M&^ \ 

8/3 4 [2ifr^sin(^)] 2 [h+2ifr^cos(V0] 2 
S 4 (r,^) 



1- ^ V^— v^j ^.^.^—^j ^ (no) 



o 7r 1 /o[S(r,^)]J I 2 [S(r,^)] 

(oob(40)>, fdr re"*' f ^ (l-^^l (j^-^f 
' ' ' tt 1 /,[H(ri)]|il 2 [S(^)] /o[H(r,V)] 

8/3 4 [2^r^gsin(i/')] 2 [/i+2A> % /gcos(V')] : 



Apparently, for rotation-invariant states we obtain the trio {Ai,A2,A3} from 
{pi, p2, P3} by making in the latter the replacement (cos(40))0 — > — (cos(40))0: 

Ai((cos(4^))^) = pi(-(cos(4^)) ) 2 = 1,2,3 (112) 

The same is true for the expressions ( |63| , |64| ) for the reflection symmetry breaking 
transition lines, which can be obtained by making the substitution (cos(40))0 — > 
— (cos(4(/>))0 in the two expressions (104.105) for those rotation symmetry breaking 
transitions which leave reflection symmetry intact. 
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6. Solution for specific choices for the pinning field statistics 

We now work out our results for specific choices of the pinning field distribution p((j>) 
with decreasing levels of symmetry, restricting ourselves (as throughout this paper) 
to those with p{4>) — p(—<f>). We start with the benchmark case h = 0, i.e. absent 
pinning fields. Wherever possible we have tested our theoretical predictions against 
extensive numerical simulations. Here the Langevin dynamics defined by equation 
was iterated with a stochastic Euler method with time step At = 10 -3 , for systems of 
size either N = 400 (with the advantage of better equilibration within experimentally 
accessible time-scales) or N = 800 (with the advantage of reduced finite size effects) . 
All simulation results shown are averages over 10 experiments. 



6.1. Absent pinning fields 



For h — the natural solution is the one with full rotational symmetry. Given the 
RS ansatz we are left with a single order parameter, q, and the two measures ( p6| , p7| ) 
become identical. Insertion of h = into (Q) and (^6|) immediately leads us to 

„-i^i 2 [3(r)] 



1 = o 



max 

9 



dr r e 2 ' 



-2f3K 2 (q- 



1 )2 _1 
2 ; 



(113) 



dr re-* r log/ [H(r)] - - log(27r) \ (114) 



with 5(r) 
^2 



2j3K^fqr. Expanding the right-hand side of (113) gives RHS(<7) = 
((3K) 2 q — -i^K^q 2 + 0(q 3 ). There is a second order transition at T c — K from 
a paramagnetic state (q = 0) to an ordered state (q > 0), and no evidence for first- 
order transitions. As the temperature is lowered further, q increases monotonically 
towards its maximum value at T = 0, as q = \ - T^/tt/AK + 0(T 2 ). Close to the 
critical point we can expand q in powers of r = 1 — T/T c and find q = + 0(t 2 ) 
(t l 0). All this is in perfect agreement with the results obtained earlier for Gaussian 
interactions |pl[ . 

For h — the RS susceptibility reduces to xrs = /3(h — q), according to (|96|), and 
thus obeys lim^^o XRS = \/k/AK and xrs = 1/2T for T > K. Close to the critical 
point, expansion in r = 1 — T/T c reveals that xrs = 1/2K + 0(t 2 ). xrs thus increases 
from the value ^/w/AK at T — to a cusp with value 1/2K at T = T c , followed by a 
monotonic 1/2T decay to zero in the regime T > T c (see also figure || ). 

Expression (|94|) for the AT instability of rotation invariant states simplifies 
similarly upon putting h = to 



" T ' 


2 noo ^ 






= dr re~2 r < 




K 







I 2 [E(r)} I 2 [E(r)} 



i 2 



7o[H(r)] I 2 [E(r)] 



1- 



lf[E(r)] 



i 2" 



(115) 



(with stability if the left-hand side is larger than the right-hand side). From this, in 
combination with the properties / n >o(0) = 0, we find upon inserting q = that the 
AT instability occurs at Tat = K, and thus coincides with the second order transition 
from a paramagnetic state to a q > state found earlier. In both models replica 
symmetry breaks as soon as we leave the paramagnetic region at T = T c — Tat = K- 
Using the properties I n (z) / Iq(z) = l—n 2 /2z+0(z~ 2 ) (\z\ — > oo) of the modified Bessel 
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functions [[30[ we can also study the behaviour of both sides of (115) for T — > 0. This 
reveals that the degree of RS instability diverges near T = 0: 

,. RHS f°°dr _ ir 2 

hm = / — e 2 = oo 

T^o LHS ./n r 



Finally we work out the functions ( 104 , 105 ) who's zeros mark contin uous 
transitions away from the rotation- invariant state. For h = it follows from ( 112 ) 
that the two types of symmetry breaking coincide, i.e. SfA = Sj c / T = E^n- For h — * 
all integrals over ip become trivial, and the quantities ( 106| , 107 ,108) reduce to 

I r d rre-^{l-&^ 



(H = - I dr re-~- r { 1- ^ffg^ }> (116) 



00 „ r r2 



... I. /f[E(r)]l f, 3/f[H(r)] 



«- 5J f o *~-- i'-^U'-w > (117) 



I) 



J 2 [S(r)] ] i?[S(r)] 
I [E(r)]j I*[E(r)] 



p 3 = I dr re-*' \ l-Wrl \ 71^737 (H8) 



For T >T C = K, where S(r) = 0, one simply obtains (pi,p2, Pa) = {\i |> 0). Hence 
T > K : Ei(T) = (T/Kf- 1 E n (T) = (T/K) 2 + - 

Close to the critical point, expansion in r = 1 — T/if reveals that (pi, p2, P3) = 
(i,i-2r,r)+O(r 3 / 2 ),S0 

t = 1-T/K: Ei(T) = 2r+C(ri) E„(T) = ^-2r+C(ri) 

Thus for T > K the rotation-invariant RS solution is stable, with in the case of model 

I this stability becoming marginal at T = T c = A', followed by restored stability for 
T < T c . Close to T = one finds (pi,p 2 ,P3) = {Ty/^/2K)(l, -1, 2) +C(T 2 ), so either 
£1,11 (T) does not exist, or S U i(T) = C(T 2 ). 

In figure § we show our theoretical results, together with those of homogeneously 
distributed pinning angles (to be studied next), by plotting the RS order parameter 
q and the susceptibility xrs as functions of temperature, and testing them against 
numerical simulations. We find reasonable agreement, given the CPU limitations on 
system size and equilibration times. In spite of the differences between models I and 

II at the microscopic level (notably in terms of frustration properties of spin loops), in 
the absence of pinning fields there is no macroscopic distinction between their physical 
behaviour in the replica-symmetric state, probably due to the overruling amount of 
frustration. The models have identical q ^ and RSB transition lines and identical 
values of the replica-symmetric physical observables, with a paramagnetic state for 
T > K, and a spin-glass state for T < K. The only difference is the degree of stability 
of the rotation invariant state against non-rotationally-invariant fluctuations (since 
Si(T) ^ Yj\\(T)) : which cannot be measured directly. 

6. 2. Homogeneously distributed pinning angles 

Our second choice is the homogeneous pinning ang le distribution p(^) = {lit)- 1 . The 
Hamiltonian (0) is no longer invariant under simultaneous rotation of all spins, but we 
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Figure 2. Theoretical predictions for the RS order parameter q (left) and the 
susceptibility xrs (right) as functions of temperature, for models I and II and with 
K = 1, in the case of homogeneously distributed pinning angles. Different curves 
correspond to different values of the pinning field strength h, taken from {0, 1, 2, 3} 
(lower to upper in left picture, upper to lower in right picture). Connected 
markers: results of numerical simulations with N = 400 (model I: •, model II: o). 



still expect the system globally to have full rotation symmetry. We cannot simplify 
equations (|8^]86|]95|) further: 



Q = 



dip , 
— / dr re 
2tt 



f[q] = -2/3K 2 (q-^) 2 



Xrs 



27r Jn 



dr 



(119) 

dr re~^ 2 log I [E(r,ip)} (120) 
J 2 |H(r,V0] ifHr^)}' 



7o[S(r,V)] / 2 [S(r»] 



h 2 +AK 2 r 2 q cos(2?/>) +4hKr^/qcos(ip) 
h 2 +4K 2 r 2 q+4hKr^/qcos(i(;) 



(121) 



with S(r, ip) = (3 \J '4qK 2 r 2 + h 2 + 4Khr^/qcos(ip). For h ^ one no longer expects 
to find a phase transition from a<7 = 0toaq>0 state; this is clear upon 
expanding (119) for (3 — > 0, giving q = \{j3h) 2 + C(/3 3 ). For weak fields one has 
q = h 2 /8(T 2 — K 2 ) + 0(h 3 ), in the regime T > K. Expansion for strong fields, 
h — > oo, gives the leading orders 

1^2' f ^ ~ h > Xrs ¥ (h^oo) 
For T — > 0, on the other hand, one finds 

1 „ r# f°° dr re- 1 



Q = 



T 



2tt 



h 2 + 2K 2 r 2 + 2KhrV2cos(ip) 



0(T 2 ) 



(122) 
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0.0 0.5 1.0 1.5 



T/K 

Figure 3. The phase diagram for models I and II, in the case of homogeneously 
distributed pinning angles. Solid: the AT instability, signaling a second order 
transition from a replica symmetric paramagnetic state (RS) to a locally ordered 
state without replica symmetry (RSB), and possibly to states without global 
spherical symmetry. The transition field strength diverges as T/K — » 0. 



lim XRS 



"dip 



dr re 



[h + V2Kr cosO)] 5 



[h 2 + 2K 2 r 2 + 2KhrV2cos{^)] 3 / 2 



(123) 



All expression reduce for h — > to those obtained in the previous sub-section (as they 
should). Away from T, h S {0,oo} we must resort mainly to numerical evaluation 
of our equations. In figure ^ we compare the result of this exercise with numerical 
simulations, carried out for N = 400 and h e {0, 1, 2, 3}. The agreement is satisfactory, 
apart from the h — and low temperature results, where finite size effects and 
equilibration problems are most prominent. We note, in comparison with the phase 
diagram of figure |l that the the most serious deviations, observed in the susceptibility 
curves, occur in the RSB region where xrs is indeed not expected to be correct. 



Next we turn to our expression (124) for the AT line 



For h 



= \' A ±\ dr 



J 2 [5(r,V>)] JflSfcVQ] 
Io[H(r,V)] ^o 2 [2(r,V)] 



1- 



I?[Z(rM 



(124) 



we recover (115) with replica instability for T < K. Expanding (124) 
for weak fields in the regime T > K gives (T/K) 2 = 1 - h 2 /2(T 2 -K 2 ) + 0(h7h), 
which has no solution; hence weak fields strengthen replica stability for T > K (as 
expected). For h — > oo we find H(r, ip) — > /3h, as a result of which we g et, u sing 
I n (z)/I {z) = l-n 2 /2z + 0(z- 2 ) (\z\ -> oo), for the right-hand side of jl24| ) the 
asymptotic form RHS = T 2 /2h 2 + 0(h~ 3 ) (h — > oo). Hence for every nonzero 
temperature there is a pinning field strength above which replica symmetry holds. 
Investigation of the limit T — > in ( |124| ), upon again using the asymptotic forms 
of the modified Bessel functions inside integrals, we now end up with the divergent 
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expression 



RHS 

lim - ; — 

t^o (T/K) 2 



r dip 

7T 



dr re 



o r 2 + \h 2 /K 2 + y/2rh cos ijj/K 



DxDy 



(x + h/K^/2) 2 + y- 



DxDy 



a xh/KV2 



'dr 



2 I [hr/KV2] = oo 



Thus the field strength hj^(T) for which the AT instability o ccur s diverges as 
T/K — > 0. For general fields and temperatures the integrals in ( 124 ) will have to 
be evaluated numerically. This gives rise to the phase diagram shown in figure [| 

For homogeneously distributed pinning angles, where (cos(A(j))) ^ = 0, it again 
follows from ( |112| ) that the two types of rotation symmetry breaking coincid e: for al l 
T one has SJ'j^T) = S™/j(T) = Siji(T). The constituents {pi,p 2 ,P3} of ( |l0^ , |l05| ) 
become 

/ 2 [S(r,V)]J 
7?[a(r,V)] 



dr re 2? ' / — < 1 



*(i0 

7T 



P2 



dr re 



*dijj 

IT 



P3 = 



dr 



(i 



1- 



I 2 [E(r,^)]jV / 2 [S(r,^)] 
/ 2 [3(r,^)]\ i?[S(r,^)] 



/o[S(r,V)]J / 2 [S(r,^)] 



(125) 
(126) 
(127) 



Since |/ n [z]//o[z]| < 1 it is clear that pi G [0, 3], pa G [— 1, 5], an d P3 £ [ 0, 1]. Hence, 
given existence of the square root, one immediately obtains from ( 104 105| ) the bounds 
Si(T) > {T/K) 2 - 1 and Sn(T) > (T/K) 2 - §. Thus the rotation invariant state is 
guaranteed to be stable for T >K (model I) and T > |v3 (model II). For weak fields 
we may again, in the regime T > K, expand in powers of h. This gives 

1 + rJ~^ + 0(h 3 ) 



Ei(T) 



Sn(T) 



T 2 -K 2 



0(h 3 ) 



(128) 



(129) 



In leading order, weak pinning fields increase the stability of the rotation-invariant 
state for model I, but do not effect the stability for model II. For strong pinning fields 
we expand our equations in powers of ft. -1 , which gives T,i v u(T) — (T/K) 2 + 0(h~ 2 ), 
so that rotation symmetry is now stable for any temperature. 

For arbitrary pinning field strengths and temperatures the integrals in 
(125,126 127) must be evaluated numerically. This reveals that for all non-zero 
temperatures and all field strengths Siji(T) > 0, implying (in combination with 
lacking evidence of first order transitions, and within the RS ansatz) the prediction 
that the system is always in a rotation-invariant state. Thus there are no further 
RS transitions, and the phase diagram is given by figure ^. As with absent pinning 
angles, models I and II differ only in the degree of stability of the rotation invariant 
state against non-rotationally-invariant fluctuations. 
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T/K 



Figure 4. Theoretical predictions for the RS order parameters q cc and q s3 as 
functions of temperature, for model I with K = 1, in the case of pinning angle 
distribution p(<f>) = \5[4> — t7t] + hfi[<f> + \tt\ (where the theory predicts that 
qcc = qss)- Different curves correspond to different values of the pinning field 
strength h, taken from {f,2, 3} (lower to upper). Markers: results of numerical 
simulations with TV = 400 (qcc'- *, qss- o). Error bars in the measurements are of 
the order of 0.02; the deviations observed are mainly finite size effects. Note that 
the property q cc = qss is the fingerprint of rotation-invariant states. 



6.3. Inhomogeneous distributions: p{fa) = \8[4>— a] + ^<5[0+a] 

Finally we turn to pinning angle distributions of the form p{fa) — i8[<j)—a] + ^5[(j>+a], 
with a £ [0, 7r]. By varying the parameter a we can control this distribution to be 
either uni-modal or bi- modal. One now has (cos(20))0 = cos(2a) and (cos(40))0 = 
cos(4a). Hence, according to ( |S8| ) we can now have rotation- invariant solutions only 
when a £ {tt/4, 3tt/A}; in the latter two cases one has (cos(40))0 = — 1, so as 
far as rotation-invariant solutions and their stability properties are concerned the 
a £ {7r/4, 37r/4} models behave identically. For a {7r/4, 37r/4} there is generally 
only the overall reflection symmetry fa — > —fa to be exploited, and we therefore 
must resort mainly to evaluating our equations for reflection symmetry states (where 
q cs = Qcs - 0) dlggl), numerically. 

In view of the above we concentrated on the cases a £ {0, 7r/4, tt/2}, which cover 
bi- modal (a = 7r/4, tt/2) and uni-modal (a = 0) distributions, and models with 
rotation invariant solutions (a = tt/4) as well as cases without (a — 0, tt/2). We 
tested our theory by comparison with numerical simulations. The general effect of 
random pinning fields is to break the symmetry between models I and II, and the 
values of otherwise identical order parameter pairs such as {q C c,qss} and {Q cc , Qss}- 
The agreement between theory and simulations is generally satisfactory, except for 
the susceptibility when measured at low values of T and h, where one faces difficulties 
associated with equilibration and replica symmetry breaking. 
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A: p(0) = i^-i7r] + i^+i7r] 

For a = 7r/4 one still has rotation- invariant solutions to our RS order parameter 
equations. In rotation invariant st ates all results obtained earlier for p (<f>) = (2tt)^ 1 
relating to the ord er pa rameter q ( 119 ), t he f ree energy per oscillator (120), the RS 
susceptibility xrs (121) an d the AT line (124) can be taken over without alteration. 
The difference with p(<f>) = (2tt)^ 1 is in the stability properties against fluctuations 
which violate rotational symmetry. Since here (cos(47r))0 = —1, we now obtain the 
symmetry breaking transitions as the zeros of the following functions (using (|112|)): 



srcr) = (T/Kf- J{ Pl {-i)- P2 (-i)Y-pi{-i) Pl (-i) - P2 (-i) 



Zf(T) = (T/Kf- v /( pi (+l)-p 2 (+l))2-^(+l) - Pl (+1) -p 2 (+l) 
EJf*(T) = {T/Kf- ^(p 1 (-l)-p 2 (-l)f-pl(-l)+p 1 (-l)+p 2 (-l) 
££ f (T) = {T/Kf- ^{ Pl {+l)-p 2 {+\)Y-pl{+\) + Pl (+1) + P2 (+l) 



(130) 
(131) 
(132) 
(133) 



with 



dr 



dr 



1- 



/ l I 2 [5(r»] 
I 2 [S(r,^)] 

_ h , 2 [*drp ( 4[gMQ] I 2 2 [S(r,^)] 
r Uo[S(r,^)] I 2 [S(r,^)] 

8/3 4 [2iYr^gsin(i/')] 2 [/i+2^r^gcos(V')] 2 

1 s*mo 



(134) 



1 



1- 



3/x 2 [S(r,^)] 



<ir re 2 ' 



,, 2 ^ d^ \ h [S(r, ^)] II [S(r, ^)] \ / J 2 [5(r,VQ] _ 3/ 2 [5(r, VQ] 

J„[H(r,^)] I 2 [S(r,^] 



o tt Uo[S(r,V)] ig[H(r,^)] 

4 



dr re 2 ' 



e?r re 2 ' 



8/3 4 [2 J ftTr^gsin(7/;)] 2 [/i+2is:r % /gcos(V')] z 
/ 2 [S(r,V0n /! 2 [H(r,V)] 



(135) 



7o[S(r,^)]J / 2 [2(r,^)] 



T # f J 2 [5(r, VQ] 1 f J? [S(r, y>)] 2 J 2 [5(r, VQ] \ 
^ /o[3(r,V)]J l^[H(r,V)] Jo[3(r,V0] J 
8(3 4 [2Kr^/q sin(^)] 2 cos(V>)] 2 



(136) 



The inequality |J n [z]//o[z]| < 1 allows one to obtain the crude general bound T/if < 3 
for the zeros of ( p^0fllll]Jl32]Jl3^ ). For strong pinning fields one has 

(Pi, P2 , p 3 ) = ^(1, -1,2) + 0(h- 2 ) (h^oo) 

and thus Ej° t I ' ref (T) = (T/Kf + 0(h~ 2 ), confirming stability of the rotation invariant 
state for sufficiently strong fields. For weak fields we may again in the regime T > K 
expand in powers of h. This reveals that those terms in {pi, p 2 , pz\ which are 
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proportional to k are of sub-leading order 0(h 3 ) as h — > 0. From this it follows 
that we revert back to expressions (128, 129Q . Apparently, the leading order effect of 



switching on the pinning fields is again to stabilize the rotation invariant state. 

For intermediate values of temperatures and fields we have evaluated 



(130,131 132 133) numerically, which showed that for all non-zero temperatures and 



all field strengths one has E™^' 16 (T) > 0. This, in turn, implies (in combination with 
lacking evidence of first order transitions, and within the RS ansatz) the non-trivial 
prediction that for p((f>) = ^5\4>— T7r] + |-£[0+T7r], which is a pinning angle distribution 
without rotation invariance, the system is still always in a rotation-invariant state. This 
statement and its quantitative consequences are confirmed convincingly for model I 
by the numerical simulation data shown in ^; similar results can be shown for model 
II. We conclude that for p(<p) = \8[4> — + ^S[<j> + \ir] the phase diagram in the 
(T/K, h/K) plane for both models I and II is simply identical to that of homogeneously 
distributed pinning angles, as shown in figure ||. 

B: p(cj>) = 5[cj>] 

For p((j>) — S[4>] we no longer have rotation invariant solutions, and the analysis 
becomes more complicated. Given the assumption of overall reflection symmetry 
(Qcs = Qcs = 0, whose stability we will calculate below) we are left with three RS 
order parameters, {Q cc ,q cc ,q ss }, and with the effective measures ( ^3|j44|) which now 
become 

Mj(9\x, y) = e P hco ^ e )+{0Kf{2Q co -i+q tB -q ac ] cos(2e)+2i3K[cos(e)x^q—+ s in(e)y^q-] 

(137) 

Me(6\x, y) = e hcos ( e )-(i 3K ) 2 [ 2 Qcc- 1 +<l»-Q^} cos(28)+2l3K[ C o S (e)x^q77+ S in(6)y^q^] 

(138) 

The remaining order parameters are to be solved from the coupled equations 

Q cc = «((cos 2 (0))4) q cc = ««cos(0)>2»>* q ss = («<sin(0))2)))4l39) 

At the RS ground state we know that q cc +q ss = 1 and Q cc = ^[l + q cc — q ss }. Due to 
p(4>) — 8[4>] the RS susceptibility (p^) simplifies immediately to: 

XRS =13(1 -Qcc -qss) (140) 

For strong pinning fields one finds lim^oo Q cc = limh^oc qcc = 1 and lim/j^oo q ss = 0. 
Hence lim/j^oo xrs = 0. In the high temperature regime an expansion in powers of (3 



shows the solution of (139) to behave as 



Qcc - \ + 0((3 2 ), q cc = lfh 2 +<D(p 3 ), q 8S =0((3 3 ) (141) 

Hence xrs = 1/2T + 0{T- 3 ) as T -> oo. 

The expressions for the AT line(s) ( p3| , p4[ ) and reflection symmetry transitions 
(^3],[34|) cannot be simplified further, except in special limits and in the special case 
described below. For h — > oo one trivially extracts from these equations that both 
replica symmetry and reflection symmetry are stable for all finite temperatures; the 
same is true for high temperatures and arbitrary field strengths (as expected). 

However, analytical progress can be made for model I. Here one observes, in 



line with (141), that the order parameter equations allow for solutions with q ss = 0. 



Given the identification q ss = lim^v^oo N ( sm ($i)) 2 such solutions imply that 



Frustrated coupled oscillators with random fields - October 23rd, 2002 



26 




Figure 5. The phase diagrams for models I (left) and II (right), in the case 
of p{4>) = S[(j>]. Solid: the AT instability, signaling a second order transition 
from a replica symmetric state (RS) to a state without replica symmetry (RSB). 
For model I the AT line coincides with line marking a continuous bifurcation of 
q ss 7^ solutions of the RS order parameter equations. 



(sin(6*i)) = for each oscillator. This can be understood as a result of the action of the 
pinning fields, which for p(<f) = 8(<f) tend to drive the oscillator phases towards 4>i = 0. 
Insertion of q ss = as an ansatz into (137) gives (sin^))* = (sin(0) cos(6*))* = 0, 
confirming q ss = self-consistently (this is not possible for model II). This results in 
considerable analytical simplifications, such as 7 CS = and j ss = (sin 2 (9))+ 7 and also 
reduces the computational effort in the various numerical integrations. One now has 
only two RS order parameters, Q cc and q cc , to be solved from 

Qcc = JDx (cos 2 (#))„ q cc = JDx (cos(#)) 2 (142) 

with the simplified effective measure 

Mi(9\x) = e P c ° s ( e )l h + 2Kx V 7 i^\+{PK) 2 [2Q c:c -i-q ac }co S {26) ^ 143 -j 

The ground state has Q cc = Qcc = 1- The equation for the AT line of model I becomes 

(T/2K) 2 = max (^Jdx [(cos 2 ((9)K- (cos((9)) 2 ] 2 , J Dx (sin 2 ((9)) 2 J (144) 

Similarly we can for q ss = simplify the constituents {Xi} of ([33]) to 



Ai = 2 J Dx (sm 2 (0) cos 2 (0))+ (145) 
A 2 = 2 Jdx (sin 2 (0)) i [(cos 2 (0)) i -2(cos(0)) 2 ] (146) 
A 3 = 4 I Dx (sin 2 (60cos(6>)Mcos(60)* ( 14 ?) 
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These are to be inserted into (|53|), whose zeros mark reflection symmetry breaking. 
Finally we now have a third type of transition: de-stabilization of q ss — 0. The 



condition for this, d //[. . ]/dqi s = 0, can with help of Appendix A be written as 



(T/2Kf = jDx{s\n\9))l (148) 

Clearly, for h — > (absent pinning fields, as studied in the first part of this section) 
the q ss ^ bifurcation occurs at T — K, and coincides with a bifurcation of q cc =/= 



and with the AT line. Comparison with (144) shows that the critical temperature T ( 



defined by (148) obeys T c < Tat- Hence the AT instability does occur for q ss = 



and is thus given by (144). Numerical analysis reveals that along the line (144) one 
always has JDx [(cos 2 (6>))*-(cos((9))2] 2 < JDx (sin 2 {6)) 2 , and hence T AT = T c (the 
two transition lines coincide, in agreement with p6|, p7| ). 

Numerical solution of equations (^,^) and ( |1 48] ) lead us to the phase diagrams 
shown in figure |[ The q ss ^ transition for model I coincides with the AT line. We 
observe that the RS stabilizing effect of the pinning fields is much greater for model II 
than for model I, as reflected in a significantly lower RSB transition temperature for 
any nonzero field strength. Numerical evaluation of the conditions ( ]63| , |64| ) marking 
the continuous breaking of reflection symmetry shows for both models that this does 
not happen at any finite temperature; reflection symmetry remains locally stable. 

The results of comparing the solutions of our RS order parameter equations 
(obtained numerically) with simulations are shown in figures ^| (model I) and ^ (model 
II). It is clear from the figures ^ and ^ that, due to the breaking of global rotation 
symmetry by the external pinning fields, there is a now a non-trivial difference between 
the macroscopic behaviour of our two models. The energetic conflicts in the two 
models are resolved differently. With the exception of the susceptibility, the agreement 
between theory and numerical experiment is good, taking into account finite size effects 
which lead as usual to a smoothening of the second order phase transition marking 
q ss 7^ bifurcations. Comparison with the locations of the AT lines in figure |s| suggests 
again that the serious deviations in the susceptibility are due to replica symmetry 
breaking; \RS cannot be expected to be correct in the RSB region. 

C: p{4>) = \S[<l>-\ir\ + \S[<l>+\n] 

Our final example is a bi-modal pinning angle distribution which, due to (cos(2#))0 = 
— 1^0, again does not allow for rotation invariant states. According to transformation 
(||) (with 7/ = A = 0) this choice is related to the case p{4>) = 8[4>] by a simple gauge 
transformation, and our models must have identical free energies and phase diagrams. 
As a consistency test we will try to extract this property from the RS saddle-point 
equations. Given reflection symmetry, the two effective measures (^3|,Q) now become 

Mi(0\x, y, ±-) = e ±! 3hsia ( e ) + (PK) 2 [2Q cc -l+q ss -q cc ] cos(20)+2/3if [cos(0)xv?^+sm(6/) av /9l7] 

(149) 

Mt!_(0\x, y, ±-) = e ±P hsia (e)-(PK) 2 [2Q cc ^l+q ss ^ qi:c ] cos(2e)+2/3if[cos(e)x A /S7+sin(6l)3/ A /9^] 

(150) 

We observe that our models can indeed be mapped onto those obtained for p(4>) — 5[4>]. 
We introduce the transformation = ±(571" — 0') (permutation and reflection of 
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Figure 6. The theoretical predictions for the RS order parameters {q S s, Qcc, Qcc} 
and the RS susceptibility XKS as functions of temperature, for model I with K = 1, 
in the case of p((f>) = 5 [</>]. Thick solid lines: theoretical predictions. Markers: 
results of numerical simulations with N = 800. Different curves correspond to 
different values of the pinning field strength h, taken from {1, 2, 3} (indicated by 
squares, circles and triangles, respectively). 



axes), which, via the order parameter equations (13E), induces a corresponding 
transformation of the order parameters: 

Qcc = l-Qcc! 1cc = q' ss , Qss^q'cc (151) 



The effect on the measures (149.15C) is that the latter can now be expressed in terms 
of the measures ( p. 37 , [138 ) describing the distribution p{4>) = 8[(f>] studied previously: 



M l (9\x,y,±-) = M I (9'\±y,x) 



M n (9\x,y,±-) = M II (9'\±y,x) 



Since operations such as (x, y) — > (±y, x) have no physical consequences, and since 
all observables and transition lines are (within RS) constructed from the measures 
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Figure 7. The theoretical predictions for the RS order parameters {q S s, Qcc, Qcc} 
and the RS susceptibility XRS as functions of temperature, for model II with 
K = 1, in the case of p(<f>) = S[<f>]. Thick solid lines: theoretical predictions. 
Markers: results of numerical simulations with N = 800. Different curves 
correspond to different values of the pinning field strength h, taken from {1, 2, 3} 
(indicated by squares, circles and triangles, respectively). 



(149,150), we conclude that the macroscopic physics of the cases p(<f>) 

J. 1 1 i liu. 



5[4>] and 

p{4>) = ^6[4>— \tt] + !<5[0+^7r] are indeed identical. The relations (151) map the 
order parameters {Q C c, Qcc, Qss} of the present bi-modal distribution to the parameters 
{Q'cc q' cc , q' ss } of the case p(4>) = S[4>], and both cases have identical free energies for 
all values of {T/K,h/K} and identical phase diagrams. In particular, the q ss = 
solution of model I for p(cj>) — 5[cj>] corresponds here to a solution with q cc = (again 
for model I only), and the RS susceptibility now becomes xrs = /?(Q CC — q cc )- The 
equivalence of the p{<j>) — S[4>] and p{<j>) — ^S[4>— iir] + ^S[(f>+^Tr] models is also 
immediately obvious in numerical simulations (which we do not show for brevity) . 
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7. Discussion 

We have studied models of frustrated coupled Kuramoto oscillators in the presence 
of random pinning fields. To simplify our problem we have assumed that the natural 
frequency of all oscillators are the same and that the pinning fields distribution is 
reflection-symmetric, i.e. p(cj>) — p(—<j>) (a simple gauge transform allows one to map 
the present models also to those with larger families of pinning field distributions). 
We have calculated the disorder-averaged free energy using the replica method for two 
types of random pair interactions (model I and model II, the first we standard ferro- 
and anti- ferromagnetic interactions, the second with chiral interactions), which differ 
in the level of frustration. In terms of e.g. oscillator triplets, model I is partially bond- 
frustrated, whereas model II is fully bond-frustrated. For small system sizes the two 
models have significantly different ground state properties. Our main interest here 
is to understand the difference(s) in physical behaviour between the two models in 
equilibrium, in the thermodynamic limit. At the mathematical level our two models 
differ in the details of the effective (replicated) single-spin measure, but they have 
the same types of order parameters. Within the replica symmetric ansatz, one finds 
different expressions for the AT instabilities; more unexpectedly, the two models are 
also found to differ in the nature of the extremum of the free energy. 

We inspected the effects of different types of symmetries which our models could 
inherit from the pinning field distribution. Global reflection symmetry leads to a 
simplified measure and a reduced number of order parameters. Reflection symmetry 
breaking transitions are found to be possible at most in the RSB regime. In the 
case of global rotation symmetry, only one order parameter remains (within RS), and 
the effective measure, the free enrgy, and the AT lines are identical for both models. 
We analyzed our order parameter equations in detail for four specific choices for the 
pinning fields, and tested our theoretical predictions against computer simulations. 
Without pinning fields, our two models both behave identically to a standard long- 
range XY model with non-chiral Gaussian interactions |3l| . The same is found to 
be true for uniformly distributed pinning angles. The agreement between theory and 
simulation is quite satisfactory, modulo finite size effects, except for the susceptibility 
at low temperatures, where RSB effects in combination with equilibration difficulties 
play a major role. We then studied inhomogeneous pinning angle distributions, of the 
form p((f>) = ^8[4>— a] + hS[(j)+a]. For the special value a = J one finds, remarkably, 
that our models are again always in a rotation invariant state (in spite of the fact that 
the pinning field distribution does not have rotation invariance); this is confirmed by 
simulations. For a — the symmetry between our two models is broken. They now 
differ significantly in terms of the location of the AT line, which for model I in addition 
coincides with a further order parameter bifurcation, as in j?6| [27j| . The a07r2 case 
can be mapped onto the a = one, by a suitable gauge transformation. Again, there 
is a good agreement between theory and simulations. 

As a next stage it would be interesting to explore the behavior of these models for 
oscillators with random distributions of natural frequencies, as well as in the presence 
of more general chiral interactions. In both extensions/generalizations the standard 
equilibrium replica formalism can no longer be used; in the first case one has to 
rely on dynamical formalisms |^] , whereas in the second case one has to call upon 
renormalization tricks IpOjl . 
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Appendix A. Derivatives of the RS free energy 



The (first and second order) derivatives of /[. . .] in (f30|) with respect to the RS order 
parameters are calculated by working out the following general relations (where 7 and 
7' denote any two order parameters from the set {Q**, 9**}): 

1 df dU 1 / II d\ogM{. ..), 



(3K 2 <9 7 <9 7 (3 2 K 2 \ \- <9 7 
1 d 2 f d 2 U 1 / II .d 2 log M{...) 



j3K 2 d^di djdy f3 2 K 2 \V d ~f d Y 

<91ogM(...)dlogM(. ..) d log M (...) d\ogM(. 



dj 97' 97 97' 

The Gaussian disorder variables {x,y,u,v} generated by the differentiations are 
eliminated via integration by parts. Below we give the final results of the calculations, 
which, although not fundamentally complicated, can be lengthy. 

Appendix A.l. First order derivatives 

For model I one finds the following first order derivatives: 
1 dh 



f3K 2 8Q CC 
1 dh 

I3K 2 dQ cs 
1 dh 



!3K 2 dq cc 

(3K 2 dq ss 
1 dh 



4{Q cc -((((cos 2 (#))4) } 
4{Q cs -<«(sm(0) cos(0))*»M 
2{(({{cos(6))l))U-q cc } 
2{{({{fw(e))l)))4- qaB } 



I3K 2 dq cs 
For model II one finds 



4{((«cos(0)Msm(0))4>0-fe} 



' <>Jn --4{Q cc -((((cos 2 (fl)M)) } 



PK 2 8Q CC 
1 dh 

I3K 2 0Q C( 
1 dh 

(3K 2 dq cc 

1 dh 

PK 2 dq ss 

1 dh 



-4{Q cs -««sin(0) cos(0))*)»4 
2{((((sm(0))l)))^q S s} 
2{{{{{co S (d))l))U- qcc } 
-4{((((cos(0)),(sin(0)),))) - gcs } 



/?^ 2 9 9cs 

The above results re-confirm our replica-symmetric saddle-point equations. 
Appendix A. 2. Second order derivatives for model I 

I^qqI = 4 - 4O3if) 2 ((((cos 2 (20)K - <cos(2tf)} 2 ») 



Frustrated coupled oscillators with random fields - October 23rd, 2002 33 
; i>2jl =4-4(/3if) 2 ((((sin 2 (20)K-(sin(2^))) 



m 2 dQls 

1 d 2 /i 
f3K 2 dQ cc dQ c 

1 9 2 /i 
f3K 2 dQ cc dq cc 

1 d 2 /i 

1 9 2 /i 

(3K 2 dQ cc dq cs 

1 9 2 /i 

/3X 2 dQ cs dq cc 

1 9 2 /i 
1 d 2 /i 

1 d 2 h 



f3K 2 8q 2 cc 




1 a 2 /: 




(IK* Oql 




1 a 2 /: 




/3X 2 Oql 





= - 4(/?if) 2 ((((sin(20)cos(20))* - (sin(20))*<cos(20)>*>»* 

= 8(/3X) 2 ((((cos(20)co S (0)).(cos(0)>, - (cos(20)K(cos(#)) 2 ))) 

= 8(/3X) 2 ((((cos(2e)sin(0)),(sin(0)), - (cos(20))*(sin(0)) 2 ))} 

= 8( / 3X) 2 ((((cos(20)cos(0)).(sin(0)>, + (cos(20) sin(0))*(cos(0))* 

-2<cos(20))*<sin(0))*(cos(0))*})}^ 
= 8(/3X) 2 ((((sin(20) cos(0))*(cos(0))* - (sin(20))*(cos(0)) 2 })) 

= 8(/3X) 2 (((( S in(20)sin(0)>,(sin(0)), - (sin(20)>,(sin(0)) 2 ))) 

= 8(/3if) 2 ((((sin(2(9) cos(0))*(sin(6>))* + (sin(20) sin(6>)}*(cos(<9)}* 
-2<sin(20))*(sin(0))*<cos(0))*})) 

- 2 + 8(/?tf ) 2 ((([(cos 2 (tf)K - (cos(0)) 2 ][(cos 2 (0)), - 3(cos(fl)} 2 ]))) 

- 2 + 8(/3i0 2 («[<sin 2 (fl)K - <sin(0)) 2 ][<sin 2 (fl)K - 3<sin(tf)) 2 ]») 



5. 

1 d 2 h 
(3K 2 dq cc dq 

C. 

i d 2 h 

(5K 2 dq ss dq c 



-4 + 8(/3X) 2 (((Ksin 2 (0)),-(sin(0)) 2 ][(cos 2 (0)),-3(cos(0)) 2 ] 
+ [(co S 2 (0)),-(cos(0)) 2 ][( S in 2 (0)),-3( S in(0)) 2 ] 
+ 2[<sin(0) cos(0))*-<sin(0))*(cos(0))*][<sin(0) cos(#))*-3(sin(#)}*(cos(0))*]))) 

- 8(/3X) 2 ((([(sin(0)co S (0)),-(sin(0)),(cos(0)).] 

x [(sin(0)cos(0)).-3(sin(0)),(cos(0)),]))) 
8(/3K) 2 ((([(sin(0)co S (0)),-( S in(0)),(cos(0)).][(co S 2 (0)>,-3(co S (0)) 2 ] 

+ [(sin(0)co S (0)).-3(sin(0)),(cos(0)),][(cos 2 (0)).-(cos(0)) 2 ]))) 
8(/3K) 2 ((([(sin(0) cos(0)),-(sin(0)),(cos(0)),][(sin 2 (0)),-3(sin(0)) 2 ] 

+ [(sin(0)co S (0)).-3(sin(0)),(cos(0)),][(sin 2 (0)),-(sin(0)) 2 ]))) 
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Appendix A. 3. Second order derivatives for model II 

' ° 2 Ui = -4-4(/3X) 2 ((((cos 2 (20)),-(cos(20))2))) 



fiK 2 dQl c 



' 1)1 h = -4-4(/3X) 2 ((((sin 2 (20)),-(sin(20)) 2 ))), 



(5K 2 dQl 

1 d 2 /n 
(3K 2 dQ cc dQ c 

1 d 2 h 

f3K 2 dQ cc dq cc 
1 d 2 h 

PR* dQ cc dq ss 
1 d 2 h 

f3K 2 dQ cc dq cs 

i d 2 h 

f3K 2 dQ cs dq cc 

1 d\h 

(5K 2 dQ cs dq ss 

1 d 2 h 

(iK 2 dQ cs dq cs 

1 d 2 h 



: - 4(/?/f) 2 ((((sin(20) cos(20))* - (sin(20))*<cos(20)>*>»* 

: - 8(/3if) 2 ((((cos(20)sin(0)),(sin(0)), - (cos(20)),(sin(0)) 2 ))) 

: - 8(/3if) 2 ((((co S (20)cos(0)),(cos(0)), - (cos(20)),(cos(0)) 2 ») 

: 8(/3X) 2 ((((cos(20)cos(0)).(sin(0)). + (cos(2#) sin(0))*(cos(0))* 

-2(cos(20))*<sin(0))*<cos(0))4)^ 
: -8(/3if) 2 ((((sin(20) S in(0)),( S in(0)), - (sin(20)),(sin(0)) 2 ))) 

: - 8(/9AT ) 2 ((((sin(20) cos(0))*<cos(0)>* - (sin(20))*(cos(#)} 2 ))) 

: 8(/3X) 2 ((((sin(20) cos(0))*<sin(0))* + (sin(20) sin(0))*(cos(0))* 
-2<sin(20))*(sin(0))*<cos(0))*))} 

: 8(/3X) 2 ((([(sin 2 (tf)K - (sin(0)) 2 ][(sin 2 (tf)K - 3(sin(tf)) 2 ]») 

: 8(/3K) 2 ((([(cos 2 (0)), - (cos(fl)} 2 ][(cos 2 (^ - 3<cos(fl)) 2 ]))) 

:4 + 8(/3if) 2 ((([( S in 2 (0)),-(sin(0)) 2 ][(co S 2 (0)),-3(co S (0)) 2 ] 

+ [(co S 2 (0)>,-(co S (0)) 2 ][( S in 2 (0)).-3( S in(0)) 2 ] 
+ 2[(sin(0)cos(0)>.-(sm(0)),(cos(0)),][(siii(0)co S (0))^3( S in(0)).(co S (0)).]))) 

1 ' h/u : -2 + 8(/3X) 2 ((([(sin(0)cos(0)),-( S in(0)),(cos(0)),] 

x [(sin(0)cos(0)).-3(sin(0)),(cos(0)),]))) 
: -8(/3/f) 2 ((([( S in(0)co S (0)),-(sin(0)).(cos(0)>.][(sin 2 (0)),-3(sin(0)) 2 ] 
+ [(sin(0)co S (0)),-3(sin(0)),(co S (0)),][(sin 2 (0)),-(sin(0)) 2 ]))) 
: -8(/3^) 2 ((([(sin(0)cos(0)),-(sin(0)).(cos(0)>.][(cos 2 (0)),-3(cos(0)) 2 ] 
+ [(sin(0)co S (0)),-3(sin(0)),(co S (0)),][(cos 2 (0)).-(cos(0)) 2 ]))) 



/3if 2 


dice 


1 


d 2 h 


/3X 2 


dq 2 ss 


1 


d 2 h 



fiK 2 dql 



(3K 2 dq cc dq 

1 d 2 h 
(3K 2 dq cc dq 

C. 

1 d 2 h 

(3K 2 dq ss dq c 
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Appendix B.l. Calculation of the RS Hessian for replicon fluctuations 

We calculate the Hessian of the disorder-averaged free energy per oscillator by 
expanding (^) in powers of the fluctuations {6q**, 8q**} around the replica-symmetric 
saddle-point. In this appendix we will use pairs of Roman indices (ab) and (de) to label 
the four combinations {cc, ss, cs, sc}, and abbreviate the corresponding functions in 
the obvious way as a(0) ,b(9) , c(0) , d(9) € {cos(0), sin(0)}. We put <j** = 2i(BKfk^ 
and 8f[. ..]=/[...] — /rs[. • ■], and obtain (since linear terms are absent): 

" w* 57[ - ■ ] = 2sin2 ( A *) E + cos 2 (^) J2 E^d 2 

af3 ab a/3 

~ af, qaf> WF\ g l fde M(e\{k**}) i/ + -'- (R1) 

with 

M(9\{k**}) = ^E e o » | ( , "^)+WE <rt E^*ft a ( , «) l ( ( f) 

Working out the fraction in the last term of ( |B.l| ) gives 

Jjg M(g|{fc**+Ms**}) 9r ^ ftn 2yv » a6 J^^(g|{fc**})°(gq)^) 

/dfl M(0|{fc**}) t^t^ / d0 M(0|{fc**}) 

. 9f/9Rn 4 v c, a6 » de /de A/(e|{fc**})«(gq)^Mg 7 )e(g g ) , 

a/3j6ab,de J y n JJ 

Thus, upon expanding log(x) = x — ^x 2 + . . ., and upon introducing the short-hand 
M(0|RS) = M (0| {fc R g}) , we arrive at 

~W SJl ' ' ] = 2s[n2{A * ] E [^Xp-5q s a c p5q c a s p\ + COS 2 (A*) £ 

afc a/3 a/3 7 5a&,de \ J M (tf 

JdB M(0\RS)a(d a )b(6 ) Jd9 M(9\RS)d(9 7 )e(0 s ) 
JdB A/(0|RS) JdO A/(0|RS) 

We next work out the last of the above quadratic terms for fluctuations around the 
RS solution JlS|), for small n. We use the short-hand Dxyuv = DxDyDuDv, and 



denote by M(0) either the measure (f26|) (for model I) or (27) (for model II): 

/ JPxyuv[l\ x Jd0M(0 x )]a(0 a )b(0 p)c(0^d(0 s ) ' 
~ \ jDxyuv[Jd0M(9)Y 

J Dxyuv[U x J d0 x M(0 x )}a(0 a )b(9 f3 ) J Dxyuv[Y[ x J d0 x M(0 x )}c(0^)d(0 s ) 
J Dxyuv[J d0M(0)] n J Dxyuv[J d0M(0)] n 

If we restrict ourselves to replicon fluctuations, where Sq^* a = Sk** a = 0, J2 a 3k a = 
J^a^QcP = an( ^ T^p^kap = J^pdQap = 0, we can proceed by inserting a string 
of Kronecker symbols (and complementary symbols 5 a p = l — 8 a p) to streamline the 
bookkeeping of possibly identical combinations of replica indices: 

SasSfjj + 8 aj dps + SasSfJj + <5q 7 <5/3<5 + 8 a ^SfjS + SasSp-/ + 5 
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giving, for n — > 0, and upon using again the replicon properties which allow us to drop 
contributions which will not survive the summation over (a, /3; ab) and (7, 5; de): 

= 8 aS 8 fil {{{[{a{6)e{6))*-{a{6))M0))A [<&(#)«*- <&(0)>*<d(0)>*]»>* 

+ Sps&cn ((({(b(6)e(6)U-(b(9)Ue(9)U [< a (0)d(0)>,- (a(0))*<d(0)>*])»* 

These two remaining terms will give identical contributions to the fluctuations around 
the RS free energy, after the summation over (7, 5; de) has been carried out, and the 
result can be written in the compact form 

- ■■] = EE \ S C S C [sin 2 (A*)C aMe + cos 2 (A*)W] 



- 25q aP 5k% 5 ab , de + \(pKY8k%Sk%E ahAe \ + . . . (B.2) 
with the two 4x4 matrices (note: a&, de G {cc, ss, cs, sc}) 

&ab,de ^ab.ss^de.cc ^ab,cc^de,ss ^ab,cs^de,sc ^ab, sendees (^**^) 

E abtde = {({{(b(e)e(9))*-(b(6))4e(9))*} [<a(0)d(0)>*- (a(0))*(d(O)) *])))* (B.4) 



Expression ( |B.2| ) shows that the replica indices (a, (3) have become irrelevant labels, 
and that for n — ► and within the sub-space of replicon fluctuations {^fc*^, the 
spectrum of the Hessian reduces to that of the following 8x8 matrix: 

(cos 2 (A*)l+ sin 2 (A*)C -I \ 

(apart from an overall multiplicative constant), with the building blocks C = {C a b,de} 
and E — {E abt de} as given in (B.3,B.4) and with the 4x4 unit matrix Iafc.de = $ab,de- 



Appendix B.2. Replicon instabilities 

Requiring the Hessian to have a zero eigenvalue (replicon instability), using the above 
results, leads us to the following condition: 

fx\ f0\ ( cos 2 [A*)l+ sin 2 (A*)C -I \ f x\ _ f 

Uy W' V - 1 4({3K)2E ) [y ) - [O 

fx\ f0\ f [cos 2 (A*)l + sin 2 (A*)C]x = y 

Uy V°y' I x = KPKfEy 
Equivalently: 

Det {4(/3X) 2 [cos 2 (A*)l+ sin 2 (A*)C]£; - 1} = 
For our two models I (where A* — 0) and II (where A* ~ ^7r) this translates into 

Model I : Det [E - (T/2K) 2 !] = (B.6) 
Model II : Det [E - (T/2K) 2 C] = (B.7) 
where we have used the property C 2 = I. 
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In the representation where the various entries of the matrices are ordered as 
{cc, ss, cs, sc}, and with assistance of the abbreviations ( |2S| . |2g| ), we find our matrices 
C and E (which will generally not commute) to acquire the following form: 



1 







-1 

0/ 



/ 7c 2 c 
7 2 

/ CS 

7cc7c 
V 7cc7c 



7 2 

/ SS 

Isslc 
Isslc 



iccic 
Isslc 
Iccls 
7 2 

ICS 



Isslc: 
7 2 

ICS 

Icclss J 



(B. 



In the representation (B^) one observes that the vector |0) = i\/2(0, 0, 1,— 1) is an 
eigenvector of both C (with eigenvalue 1) and E (with eigenvalue (((icdss — l1 s l))<t>)^ 
and therefore also of the two matrices in (B.6 B.7). This immediately leads us to the 
first replicon instability condition: 

(T/2^) 2 -(« 7cc7ss -7 c 2 s ))) (B.9) 

(for both models I and II). 

The remaining three eigenvalues of the relevant matrices in (B.6. B.7), viz. 
E— (T/2K) 2 1 (for model I) and E-(T/2K) 2 C (for model II), must all be orthogonal 
to |0), and are thus in the sub-space spanned by the following three orthogonal 
eigenvectors of C (with eigenvalues {1,-1,-1}, respectively): 



ID 



~V2(1, 1,0,0), 



|2) = 5^2(1,-1,0,0), |3) 



^(0,0,1,1) 



After some simple algebra one finds that on the basis {|1), |2), |3)} the matrices (B 
reduce to the following 3x3 ones 




(((|(7 2 C +7 2 J+7 2 S »)0 (((Ul 2 cc-l 2 s 



(((icsilc 



E 



{ih{l 2 cc-l 2 ss))))* (((hhl + lD-ll)))^ («7c(7cc-7~)>»* 

(((icsilcc + lss))))^ (((lcs{lcc-lss))))4, {{{lcclss+ll s )))<t> 



(B.10) 



(B.ll) 



The remaining replicon instabilities now follow upon inserting (B. 10, B.ll) into the 



conditions (B.6, B.7). The physical RSB transition associated with the combined 



replicon instabilities is the one occurring at the highest temperature. 



